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ABSTRACT 

We investigate the influence of large-scale stellar feedback on the formation of molecular clouds in 
the Large Magellanic Cloud (LMC). Examining the relationship between Hi and 12 CO(J=1-0) in 
supergiant shells (SGSs), we find that the molecular fraction in the total volume occupied by SGSs is 
not enhanced with respect to the rest of the LMC disk. However, the majority of objects (~ 70% by 
mass) are more molecular than their local surroundings, implying that the presence of a supergiant 
shell does on average have a positive effect on the molecular gas fraction. Averaged over the full SGS 
sample, our results suggest that ~ 12-25% of the molecular mass in supergiant shell systems was 
formed as a direct result of the stellar feedback that created the shells. This corresponds to ~ 4-11% 
of the total molecular mass of the galaxy. These figures are an approximate lower limit to the total 
contribution of stellar feedback to molecular cloud formation in the LMC, and constitute one of the 
first quantitative measurements of feedback-triggered molecular cloud formation in a galactic system. 

Keywords: galaxies: ISM, ISM: bubbles, ISM: evolution, ISM: molecules, Magellanic Clouds 



1. INTRODUCTION 

The formation of cold, dense molecular clouds from the 
atomic interstellar medium is a key process in the evolu- 
tion of galaxies, and one that sets fundamental bound- 
aries on star formation rates. While recent work has 
shown that the presence of molec ules is not strictly nec - 
essary for ongoing star formation ( Glover k Clarkj 2012), 
molecular tracers nevertheless provide the best obser- 
vational diagnostic of the coldest, densest phase of the 
ISM in all but the most metal-poor systems. Astrophys- 
ical drivers of dense (molecular) cloud formation i n disk 
galaxies include global gravitational instabilities dWada 
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matter m spiral shocks (Kim fe Ostriker| |2006[ |Dobbs 
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and compression in large-scale expan ding s hells driven 
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by stellar feedback ( Hartmann e t al. 
et al.||2002| |NtormousTet al.||2011| );m an" 



2001| |Elmegreen| 
ill cases aided by 



turbulent co mpression that enhances density on a, rang e 



of scales (e.g. |Elmegreen|2002| [Glover k Mac Low|200"7 ). 

However, disentangling the relative contributions of these 
processes is usually not trivial, and the primary drivers 
- as well as the details of the physics - are still open to 
debate. 

The role of stellar feedback in molecular gas produc- 
tion is of particular interest, since it is key element of 
how star formation self-regulates. The cumulative en- 
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ergy input from stellar winds, supernovae and ionizing 
radiation is able to form molecular clouds via the accu- 
mulation, compression and cooling of the ambient ISM 
in giant (100^1000 pc) su pershells around OB clusters 
( McCray k Kafatos |l987|). From a theoretical perspec 
tive, gravit ationaf""[e.g. Elmegreen et al' 1 "'""" ■•' 



et al. 



2012) 



1989 



flu id dynamical (e.g. Mac Low et al 
Visfm iacl |l994), and thermal (e.g. Inoue &: Inutsuka 



2002; Wiinsch 



instabilities help to concentrate and fragment the 
walls of expanding shells into dense clouds; and more 
broadly, supershells fall within the scope of systems cov- 
ered by the colliding flows theory of molecular cloud for- 
mation, in which molecular gas is formed rapidly from 
the atomic medium in the interfaces of turbulent ISM 



flows (e.g. |Audit k Hennebelle||2005| |Vazquez-Semadeni 
et al.|2006[ \'M7\ |rlennebelle etal|2O08l [Heitsch k Hart ' 
maim 2008 al|Banerjee et al.|2009||Inoue fe lnutsuka|2 012 ; 
" ark et al.||2012p . Parameter space considerations sug- 



gest that small, thermally-driven condensations will de 
velop before shells are able to become unstable on the 
longer wavelengths assoc iated with gravitational insta- 
bility (Heitsch et al. 2008 ). However, global gr avitational 



contraction helps to reduce overall timescales ( Heitsch k 



Hartmann 2008b ) , and results suggest that the onset of 
self-gravity and the conversion to the molecular phase 



likely occur at roughly the same tim e (Hartmann et al. 
2001| |Vazquez-Semadeni et al.||2007 ) 



Despite theoretical advances, there is as yet no con 
sensus on contribution of supershells to molecular cloud 
formation rates, and there is a strong need for data 
to constrain theory. Observationally, the association of 
molecular clouds with s upershells has been reported in a 
number of studies (e.g. |Jung et al.||1996| |Koo k Hciles 
19881 |Carpenter et al j[2ljOO| |Kim k Koo||200uj [McClure 
Griffit hs et al. 1 2000; Matsunaga et al. 2001; Yam agucnT 
et al.|2001a|b| [Dawson et al.|2008||2011b| ), and has often 
been interpreted as evidence of tlie m-siti 



interpreted^ 
molecular gas in shell walls (e.g. 
|fc Koo| [20001 |Matsunaga et al.| 



situ formation of 
Fukui et al.|1999aHKim 



2 001 1 |Yamaguchi et al. 
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20011)). However, while the presence of molecular clouds 
in shell walls is indeed suggestive, it is important to make 
the distinction between association and formation; some- 
thing that generally proves challenging for individual ob- 
jects. Furthermore, large-scale stellar feedback can also 
be destructive - both in the initial io nization and dy- 
nami c al disruption of birth clouds (e . g. Dale & Bonnell 
" Dale et all 120121 iWalch et all 120121), and in the 
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destruction of pre-existing molecula r clouds in the 



passage of an expanding shell (e .g. Klein et al.|1994 Pit- 
tard||2011||Dawson et aL||2011a | 
From a galactic evolution pe 



perspective, the net effect 
of large-scale stellar feedback on the molecular gas frac- 
tion is therefore of critical importance - the question of 
whether a supershell ultimately causes an increase or a 
decrease in the quantity of molecular gas in the volume it 
occupies. On e of the first studies to directly address this 
question was Dawson et al. (2011b), who found evidence 
of a net increase in the molecular gas fraction in two 
Galactic supershells. However, while the close proximity 
of Milky Way objects affords some advantages, Galac- 
tic studies suffer from difficulties in reliably identifying 
supershells and their associated emission in a confused 
Galactic Plane. To explore the question of the global 
influence of large-scale stellar feedback on a galactic sys- 
tem, we now turn away from the Milky Way, and towards 
our nearest star-forming neighbor, the Large Magellanic 
Cloud (LMC). 

In this paper we compare the molecular and atomic 
content of the LMC disk within and outside of super- 
giant shell (SGS) systems in order to investigate the in- 
fluence of large-scale stellar feedback on the molecular 
gas fraction of the ISM. SGSs are the largest and most 
energetic shells in the LMC, and those which we might 
expect to have the most dramatic effect on the evolution 
of the interstellar medium. The LMC is an attractive 
target for this kind of study, due to its relative proxim- 
ity (d ~ 50 kpc), nearly fa ce-on orientation (i ~ 35°; 



van der Marel fc Cioni|200T ) , and the presence of a large 
population of sup ershells and shell-like structures i n its 
gaseous disk (e.g . |Meaburn||1980[ |Kim et al 



1999 Chu 



& Mac Low 19901. Excellent spectral-lme data exist at 



aremmute resolutions for both the atomic and molecu- 
lar c omponents of the ISM; the former in the H i 21 cm 
line ( [Kim et al.|1998||Staveley-Smith et al.|2003[) and the 
latter in rz CO( J=l-0)lFukui et al.|1999b[ |Mizuno et al 



2001[ |Fukui et al.|2 008[). In section |2| we describe these 
datasets. Section 15] then outlines our method of select- 
ing and defining supergiant shells, describes the basic 
analysis and its assumptions, and examines sources of 
uncertainty. Section [4] describes our results, which are 
discussed further in section [5] We finally summarize our 
conclusions in section |6] 

2. DATASETS 
2.1. Hi 

The LMC was observed in the Hi 21 cm line with the 
Australia Telescope Compact Array (ATCA) and the 
Parkes 64 m telescope, both of which are operated by 
CSIRO Astronomy and Space Science (CASS). Details 
of the observations an d da ta combination are given by 
Kim et al.[ ( fl998[|2003j ) and |Staveley-Smith et al| [2003t I . 
The combined datacube covers an area of 7.5" x 7.5" at 



an effective resolution of 1' (~ 15 pc), and a la bright- 
ness temperature sensitivity of 2.4 K in a 1.65 km s _1 
velocity channel. For consistency with the CO data, the 
cube was regridded to the kinematic local standard of 
rest (LSR) velocity frame and smoothed to a resolution 
of 2.6'. The main Hi integrated intensity (/hi) image 
used in this work is integrated over the velocity range 
175 < Vlsr < 375 km s" 1 , which contains all emission 
in the disk. 

2.2. 12 CO(J=1-0) 

An unbiassed survey of the LMC in the 115 GHz 
12 CO(J=1-0) line was carried out with the 4 m NAN- 
TEN telescope in Las Campanas, Chile, operated by 
Nagoya University. The survey covered a total area of 
~ 30° and consisted of ~ 26,900 positions observed on 
a 2' grid, with an integration time of ~ 3 minutes per 
pointing. This spacing slightly undersamples the main 
beam, which has HPBW ~ 2.6' (~ 40 pc). The mean la 
brightness temperature sensitivity is ~ 0.06 K in a 0.65 
km s _1 velocity channel, but varies between ~ 0.04 and 
~ 0.09 K in different sub-regions of the map. Further 
details o f the observing stra tegy and data reduction are 
given by |Fukui et al |(|20M | 

Following I Wong et al.l (1201 



]2Dll ), we apply a smooth-and- 



mask technique to extract real emission features from 
the data. This method is is chosen because it is effective 
in recovering weak emission features from the CO cube, 
but results in no final degradation in resolution. The 
datacube is first convolved in the spatial domain with 
a Gaussian of twice the width of the telescope beam, 
smoothing it to a resolution of 6'. This suppresses noise 
which should be uncorrelated between beam areas, but 
picks out any signal which is spatially resolved by the 
beam. A mask is then generated from this cube which re- 
tains only voxels detected above the 3<r noise level. This 
mask is applied to the original unsmoothed cube, and 
the resulting emission integrated over the full cube ve- 
locity range of 200 < vlsr < 305 km s" 1 to produce 
a map of integrated intensity, Iqo- This velocity range 
differs slightly from that of Hi, but contains all detected 
CO emission in the LMC ( |Mizuno et al.|2001[ ). The mean 
RMS noise of 0.06 K per 0.65 km s _i channel corresponds 
to an integrated intensity sensitivity limit of 0.07 K km 
s _1 for a line detection across three contiguous channels. 
We perform a final cleaning step to suppress noise in 
the Jco m ap by blanking pixels below Aa = 0.28 K km 
s _1 . This results in a final map in which weak spectral 
features are well-recovered with no loss of resolution. 

3. METHODOLOGY 

3.1. Object Selection 

The disk of the LMC is rich in bubbles, shells, arcs 
and filaments, and numerous references to shells and 
supershells are present in the literat ure, identified in a 
variety of observational tracers (e.g. |Davies et al.||1976 
~^ |T980| jChu fc Mac Low||l990| |Kim eFaT|[TW 
KDSB99). However, in general there is no one- 
to-one correspondence between different studies, with 
different tracers and methods favoring different classes 
of objects. A study on LMC shells must therefore give 
careful thought to its choice of sample. 
The most comprehensive census of shells in the LMC 
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Figure 1. Supergiant shells and shell complexes overlaid on integrated intensity map of the LMC. The greyscale image is Hi integrated 
over the velocity range 175 < vlsr < 375 km s — 1 and smoothed to a resolution of 2.6'. Blue contours are 12 CO(J=1-0) integrated 
intensity, processed as described in j ]2,2| and integrated over the velocity range 200 < "UlS-B < 305 km s — 1 . (Both images contain all 
emission detected in each tracer in the LMC). Contours start at 0.65 K km s — 1 and are incremented every 2.5 K km s- 1 . The solid blue 
line marks the boundary of the region o bser ved in CO. Dark pink lines trace the inner rims of the shell complexes and purple lines mark 
their outer boundaries, as described in §3.2| Each shell complex is numbered and its constituent members listed in table [3] Dotted blue 



lines enclose the region known as the southeastern Hi overdensity (SEHO) 



disk is provided by the Hi 21 cm line, which probes 
the structure and dynamics of the atomic ISM. Super- 
shells remain detectable in Hi long after their central en- 
ergy sources have switched off, and objects defined in Hi 
should in theory comprise the most complete sample of 
shells and shell-like structures in the galaxy. The largest 
and most up-to-date catalogue of Hi supershells in the 
LMC was compiled by KDSB99. They define shells as 
coherent, bright-rimmed voids that are present in multi- 
ple velocity channels, and identify a total of 126 objects 
ranging in radius from ~ 50 to ~ 700 pc. 

In this work we consider only the 'supergiant shells' 
(SGSs) of KDSB99 - defined as those with radii greater 
than the scale heig ht of the LMC gas disk (z g ~ 180 pc). 



Book et al. ( 2008 1 compare Hi and new Ha data to ex- 
amine in detail the SGSs from KDSB99 and the 9 similar 



objects from the classic Ha catalogue of Meaburn ( 1980 1 . 
Their analysis, with which we concur, ultimately rejects 
several SGSs as false detections, and refines the positions, 
velocities and extents of some others. Our final sample 



of supergiant shells is based on this revised selection, and 
contains 19/23 of the original KDSB99 catalogue. These 
are o rganized into 11 spatially distinct complexes (see 
|3.2[ ), which between them occupy ~ 40% of the area of 
the main Hi disk. These complexes are shown in figure 
[T]and are listed, along with their constituent shells, in 
table [0 

There are a number of reasons - both scientific and 
pragmatic - for working only with supergiant shells and 
excluding smaller objects. SGSs are a physically mean- 
ingful population of the largest and most energetic shells 
in the LMC, and are therefore expected to have the most 
dramatic effect on the ISM. Their regions of influence are 
large, giving them ample opportunity to accumulate the 
large quantities of material needed for molecular cloud 
formation. They are also well-resolved at the resolution 
of the CO datacube, and far larger than a typical giant 
molecular cloud (GMC). This avoids ambiguity in the as- 
sociation of molecular gas with any given shell, and also 
avoids small-number problems that might arise from the 
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Table 1 

SGS Complexes and Thresholding Parameters 



SGS Complex Previous Listings Thresholding Parameters c Shell c (') v ex p b,g 

Ha a Hi" 6 cT 3 I 11 v LS r range c 7p Thickness km s" 1 
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LMC1 


SGS3 


04 59 00 


-65 40 00 


270-285 


15 


17 


15 


2 


LMC2 


SGS20 


05 42 45 


-69 55 00 


220-230 


35 


9 


25 




LMC2 


SGS19 


05 44 45 


-69 14 00 


260-270 


30 




25 


3 


LMC3 


SGS8/12/13 


05 28 30 


-69 17 00 


250-265 


20 


9 


23/25 




LMC3 


SGS15 


05 34 00 


-68 44 00 


265-275 


20 




24 


1 


LMC4 


SGS11/14 


05 31 45 


-66 50 00 


275-300 


15 


15 


36 


5 


LMC5 


SGS7 


05 31 45 


-66 50 00 


285-305 


20 


11 


25/30 


6 


LMC6 


SGS2 


04 59 15 


-68 40 00 


250-270 


15 


13 


17/20 






SGS5 


05 02 30 


-68 25 00 


265-285 


15 




20 


7 


LMC8 


SGS4 


05 01 00 


-70 27 00 


220-240 


10 


9 


23/15 


8 




SGS16/17/22 


05 39 00 


-68 18 00 


275-280 


20 


7 


18/26/? 


9 




SGS23 


05 51 30 


-67 37 00 


275-285 


10 


15 


23 


10 




SGS6 


05 14 30 


-65 20 00 


280-290 


10 


19 


18 


11 




SGS21 


05 45 00 


-66 30 00 


280-290 


20 


13 


19 



Mcaburnjj 


1980t 


iKim et al. 


jl999|l 



Approximate center position of thrcsholdcd void. 

Velocity range over which emission is averaged prior to thresholding. 
Threshold intensit y level. 
Book et al~1 fl2008|l 
Staveley-Smith et aT] |2003} 



association of only a few CO pixels with smaller objects. 
In addition, small/young supershells are expected to re- 
tain a correlation with the star-forming cloud complexes 
that birthed them, leading to difficulties in separating 
correlation and causality when interpreting their molec- 
ular gas fractions. SGSs, on the other hand, should be 
free of this bias, since their large size and advanced evolu- 
tionary state means that they have expanded far beyond 
their initial birth sites. In addition, SGSs are well stud- 
ied, and multi-wavelength data exists for most objects, 
meaning that there is minimal ambiguity in their iden- 
tification or origin. The same is not true of smaller Hi 
shells, many of which are likely less robust detections. 
For these reasons we restrict ourselves to exploring the 
influence of the largest and most obvious objects in the 
present work. However, w e dis cuss the implications of 
excluding smaller shells in |5.1| 

Simulations demonstrate that it is possible to form gi- 
ant dense-rimmed voids by gravity, turbulence and ther- 
mal instabilities alone, even without the inclusion of stel- 



lar feedback 


e.g. 


Wada et al. 12000 


Dib & Burkert 2005 


Dobbs et al. 


20ll 


). However, there is good evidence to 



back shells. Expansion is a strong signature of formation 
via a central energy source, and is confirmed in all of 
the constituent members of the eleven SGS complexes, 
in the form of Hi emission at velocities positive or nega- 
tive of the systemic velocity (IKim et al. |[l999; Staveley- 
Smith et al. |2003l [Book et al IjMTHj se e table |l| 



observations by Points et al.l (|200ip find diffuse X-ray 



emission in six of the eleven complexes in our sample, 
interpreted as hot gas produced by stellar winds and su- 
pernova blasts. The exceptions are Complexes 6, 7, 9, 10 
and 11. Of thes e, only Complex 6 wa s observed in their 



study; however, Bomans et al. ( 1996 ) do claim a detec- 
tion for part of this complex in the same ROSAT bands, 



Gruendl (2011 ) claim that extended gamma ray emission 
in the direction of the shell arises from cosmic ray accel- 
eration in SNR s hocks, and young clusters close to the 



edges of the shell ( |Glatt et al.|2010 ) a re also directly im - 
plicated in expansion in its JNlE lobes ( |Book et aH2 008). 
Known clusters are similarl y implicated in the dynamics 
of Complexes 1, 2, 3 and 4 (|Book et al " 



et al. 2009 Glatt et al. 2010). Of these, Complexes 2 



2008 Kawamura 



so the literature is ambiguous. For Complex 7, Chu fc 



and 4 bear special mention as objects for which a large 
body of research stretching back several decades makes a 
particularly ro bust and comprehensive case for a stellar 
origin (see e.g.J C aulet et al.|198 2] Jgppita et al.|1985||BcT 
mans et al.||199l| |Foints et aITO99p . Little information 
exists on complexes 10 and 11. However, these objects 
cover very small areas in the observed region, and their 
i mpa ct on the results of this work is minimal (see also 
{4.2). It should be noted here that Ha emission from 
shell walls is not a reliable indicator of a stellar origin, 
alt hough it is someti mes treated as such. As pointed out 
by Book et ah] ( 2008 ) , the ionizing radiation in Ha-bright 
SGSs frequently arises from young sources within shell 
walls and not from a central parent cluster. Similarly, an 
Ha-dark shell is one in which no recent (< 12 Myr) star 
formation has occurred, but not one whose origins were 
necessarily non-stellar. 

3.2. Supergiant Shell Definition 

This work compares volumes affected by supergiant 
shells with volumes that are not, using two-dimensional 
velocity-integrated intensity maps summed over the full 
velocity range of the LMC. This 2D approach is well 
justified in the case of supergiant shells. Because SGSs 
are by definition larger than the Hi scale height of the 
LMC (KDSB99), they are assumed to have propagated 
through the full vertical extent of the gaseous disk to 
blow out on both sides of the plane. An SGS is there- 
fore conceived as a cylindrical section of disk, all emis- 
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-450 -300 -150 150 300 450 

Distance from shell rim (pc) 

Figure 2. Example of the definition of shell thickness for two SGS 
Complexes. The solid and dashed lines show variation of Hi and H2 
surface density with distance from the shell inner rim, calculated by 
summing emission within 2' bands incremented inwards and out- 
wards from the thresholded inner rim of the SGS c omp lexes. Units 



are M, 



PC 



with masses derived as described in S 3.3 



The bright 



rim and evacuated void configurations can be clearly seen. Emis 
sion within the grey regions are defin ed as the SGS zones, whose 
edges is defined as described in §3.2| Uncertainties are estimated 
as ±50% of the shell thickness, and are marked with dotted lines. 
H2 is included here for illustration, but is not used in defining the 
shells. 

sion above and below of which is also related to the shell 
system. This mitigates the problems of line-of-sight con- 
tamination and shell curvature that would severely af- 
fect smaller objects. Emission at all velocities within the 
spatial extent of an SGS may therefore be assumed to 
be related to it. (The inclination of the LMC causes de- 
partures from this assumption, which we discuss below 
in 



3.4 1 



The shapes and extents of SGS regions are defined us- 
ing a thresholding method, and pixel masks are generated 
that include the entirety of the bright shells and their 
evacuated voids. This method makes no ad-hoc assump- 
tions about the morphology of a shell, instead allowing 
the configuration of the ISM to define its boundaries. 
SGSs are seen as rim-brightened voids in Hi, and are 
most clearly discerned at their systemic velocities. The 
following steps are performed to define individual SGSs 
and overlapping complexes: 

1. A narrow- velocity Hi mean intensity map is produced 
by averaging the Hi brightness over the central few ve- 
locity channels of a supergiant shell, where the exact ve- 
locity range is chosen to include the entire bright rim of 
the shell over its full spatial extent (~ 10-20 km s _1 ). 

2. A contour is plotted at a level that picks out the 
shell's inner rim. Because Hi brightnesses in the LMC 
vary with both location and velocity, there is no 'one size 
fits all' threshold level. The choice of threshold is made 



by eye, and is typically approximately half of the mean 
brightness of the shell walls. Since the analysis in this 
work is insensitive to the exact positioning of the inner 
rim, this approach is acceptable. Similarly, since many 
shells are not completely enclosed rings, it is sometimes 
necessary to 'fill in the gaps' in this rim by hand. Such 
gaps are usually small, and human judgement does not 
strongly influence the final size or shape. 
3. All pixels within this inner rim are masked, and the 
masks for all shells are then combined. At this stage, 
overlapping objects are amalgamated into single com- 
plexes. 

Because many structures designated as SGSs by 
KDSB99 are either located completely within the spa- 
tial boundaries of larger shells, or are well connected to 
them in both space and velocity, it is not necessary to 
perform the above thresholding steps for every individ- 
ual object in the KDSB99 catalogue. Where listed SGSs 
overlap significantly in both space and velocity, a single 
thresholding step at the appropriate velocity range recov- 
ers all separately catalogued features. Where there is full 
spatial overlap there is no need to threshold the smaller 
SGSs at all, since they occupy no unique pixels outside of 
their larger counterparts. Table [T] lists the approximate 
central positions, velocity ranges and threshold intensi- 
ties used to generate the final SGS complexes. 

At this stage the masks are imposed on the 2D Hi in- 
tegrated intensity map, and contours are drawn at inter- 
vals of 2' incrementing both inwards and outwards from 
the rim of each shell complex. The Hi emission between 
each pair of contours is then summed to produce a plot 
of mean surface density vs distance from the complex 
rim. Example plots for Complexes 1 and 4 are shown in 
figure [2j Low- intensity voids and high-intensity walls are 
clearly seen, demonstrating how SGSs are a dominant 
component of LMC disk structure, even when emission 
is fully integrated in velocity. The outer edge of the wall 
is defined as the distance at which the mean Hi surface 
density has fallen to the first minimum, or in cases where 
the emission falls off smoothly with no subsequent rise 
(3/11 complexes), as the distance at which the mean sur- 
face density has decreased by half of the difference be- 
tween the peak and the inner rim. It should be noted 
that the presence of a second peak outside the shell wall 
depends on the configuration of the surrounding ISM, 
and there arc no differences observed in the properties 
of the SGS systems with and without this feature. The 
mask for each shell complex is then grown by this thick- 
ness. Thus, the final set of masks cover both the bright 
shells and evacuated voids of the LMC supergiant shells. 
The locations and outlines of the final SGS complexes 
are shown in figure [I] 

3.3. The Metric: Molecular Fraction in SGS and 
'Background' Zones 

The ISM in the volumes occupied by supergiant shells 
is compared to that outside them, to seek for evidence 
of molecular cloud formation in shell walls. We aim to 
answer the question of whether the passage of an SGS 
through the ISM causes the production of more H 2 than 
would have been created if no SGS had been present. 
The underlying assumption is that if the net effect of a 
supergiant shell is to produce more molecular matter in 
its swept-up walls than is destroyed by its initial pas- 



G 



sage through the ISM, then the volumes now occupied 
by SGSs should be more molecular than they were prior 
to the occurrence of the shell. The appropriate region 
over which to measure the molecular fraction is therefore 
the entire volume affected by a supergiant shell, includ- 
ing both the over-dense walls and the under-dense void 
(which in a 2D analysis will also contain emission from 
the front and back limbs of a shell). Selecting shell rims 
alone would preferentially sample only the highest Hi col- 
umn densities, and would not probe the effect of SGSs 
on the full volumes of space that they influence. 

In the absence of information on the ISM at past 
epochs, a comparison is made between SGSs and 'back- 
ground' (non-SGS) regions that serve as a proxy for the 
undi sturbed medium. We consider two approaches. In 
§4.1| SGS volumes are compared with the rest of the LMC 
disk. In this case SGSs between them occupy 45% of the 
total observed area of the LMC, corresponding to 13.0 
deg 2 (a projected area of ~ 11 kpc 2 ), with a remain- 
ing 16.6 d eg 2 (~ 13 kpc 2 ) designated as the background 
zone. In §4. 2 1 they are compared individually with their 
local surroundings, in order to minimize the impact of 
galaxy-scale variations in the molecular fraction. In the 
latter case, background zones are defined as bands of 
constant thickness around each shell complex, excluding 
pixels that are assigned to other shells. The minimum 
background zone width considered is 10' (~ 150 pc) from 
the shell edge, and the maximum considered is 33' (~ 500 
pc) from the inner rim. 500 pc represents an approximate 
limit at which emission is considered 'local', roughly set 
by the distance at which the background zones of outer 
galaxy shells such as Complexes 4 and 8 begin to sam- 
ple the very different en viron ment of the dense central 
regions of the disk (see [4.1|. SGS complexes are ex- 
cluded from this individual analysis if either their shell 
or background zones do not include sufficient area in the 
CO datacube, defined as at least 0.25 square degrees of 
coverage (~ 0.2 kpc 2 ) in each zone. On this basis Com- 
plexes 10 and 11 (SGS6 and SGS21) are excluded. For 
this object-by-object analysis, pixels that belong to the 
swept-up wall of more than one complex (e.g. the bright 
ridge between Complexes 4 and 5) are assigned to the 
nearer of the two. Figure [3] illustrates the SGS and max- 
imum background zones for each complex. The surface 
areas occupied by SGSs and their partner background 
zones are given in table [3j 

We define two quantities, /h 2 ,SGS an d /H 3 ,bgj where 
/h 2 is the molecular fraction of the ISM, given by 



/h 2 



M fl 



HI 



Mv 



(1) 



and the subscripts 'SGS' and 'bg' indicate measurements 
in supergiant shell volumes and background regions, re- 
spectively. Mm and Mjj 2 are the atomic and molec- 
ular hydrogen masses measured from Hi and CO, and 
since the distance to all points in the LMC is assumed 
to be equal, /h 2 is simply related to column density as 
/h 2 = 27Vh 2 /(-ZVhi + 27Vh 2 ). Nh 2 is calculated assuming 
an Too ^0 Nn 2 conversion factor of 3.0 x 10 20 cm" 2 K _1 
km -1 s (appropriate for the LMC; Leroy et al. 2011 ), and 
Mhi is assumed to be equal to 1.8 x 10 iB /hi (Dickey & 
Lockman||1990[ ); the exact result for optically thin gas. 



3.4. Sources of Uncertainty 

The uncertainties inherent in estimating column den- 
sities from Hi and CO, and their potential effects on 
the observed relationship between iVm an d A^h, in the 
LMC have been discussed extensively by |Wong et al.| 
(2009 hereafter W09). For Hi the assumption of an op- 
tically thin gas will provide an underestimate of iVm by 
a factor of r/[l — e~ T ], with absorption studies towards 
background continuum sources suggesting that this fac- 
tor reaches a value of ~ 2.0 towards the most opaque 



LMC sight lines ( Dickey et al.fl994 Marx-Zimmer et al. 
20001. For CO, a, number of factors may affect the ac- 



curacy of the Nu 2 conver sion, particularly on the scales 



of individu a l clouds (e.g. [ Shetty et al. 2011a|b 



Gcnzcl 



et al 2012 |Glover fc Clark||2012[ ). These uncertainties 
are ditncuft to quantity. However, their impact on /h 2 is 
mitigated when summing over large regions of space, as 
we do here. 

The assumption that all emission along the line of sight 
to an SGS is genuinely part of the shell system implic- 
itly assu mes a face-on galaxy. For an i nclination angle of 



35° (van der Marel & Cioni 2001) and an exponen 



tial Hi disk with a scale height of 



Padoan et al. 
in a layer 



2001), 



180 pc (KDSB99; 
50% of the Hi mass is contained 



250 pc thick. A given line of sight through 
the LMC therefore contains emission that is spread over 
a projected distance of ~ 250 tani w 175 pc in the plane 
of the galaxy, translating to an in-plane positional uncer- 
tainty of ±87 pc. This will not strongly affect the shapes 
of the SGSs themselves. The thresholded inner rims iden- 



tified in {3.2 reflect a real structural property of the ISM 
- the void / rim configuration that is the signature of a 
swept-up shell. Where problems may arise is in the as- 
signing of emission to either an SGS or background zone, 
with the spreading of emission within the plane along a 
given sight line inevitably resulting in some emission be- 
ing falsely assigned. This is not an effect that can be ad- 
equately expressed as an uncertainty on the defined shell 
thickness, nor can it be well modeled even by consider- 
ing the third dimension of velocity. (Even on the steepest 
parts of the rotation curve SR = 90 pc translates to a 



Kim 



rotational velocity difference of only < 3 km s 1 ; 
etaLl[T998l |van der Marel et al.||2002l jOlsen &: Massey 
2007 ) It is a fundamental limitation of" 



) It is a fundamental limitation of an inclined sys- 
tem that cannot be corrected for without knowledge of 
the 3D configuration of the gas. Fortunately, the molec- 
ular fraction in SGSs and their background zones is most 
sensitive to the placement of the scattered concentrations 
of CO, whose positional uncertainty is expected to be less 
severe. While the vertical distribution of molecular gas 
in the LMC is not known, CO scale heights are typically 
between two and se yeral times smaller than Hi, both in 
the Milky Way (e.g [ Dickey fc Lockman| l990; M alhotra 
1994| |Stark fc Lee||2005|) and e lsewhere (e.g. |Sancisi fc 
AHen| 1979| |Scovillc et af ||1993[ ). For an exponential CO 
scale height of ~ 90 pc, the projected in-plane positional 
uncertainty on CO clouds is only ~ ±43 pc - approxi- 
mately one resolution element. It is also worth noting 
that some estimates place the characteristic thickness of 
the LMC gaseous disk even small er, at ~ 100 pc for the 
Hi layer (lElmegreen et al. 2001), which would further 



reduce this uncertainty. 

Our shell definition method relies on judging the end 
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Figure 3. SGS complexes and their local background zones, denned as described in j j3.3| Filled purp le ar eas show the SGS zones and filled 
blue areas bordered by dashed lines shows the maximum widths of the local background zones (see j]4,2|l. Spatial axes are right ascension 



and declination, and the color scale is in column density units of 10 21 cm 2 . All other lines and colors are as in Figure [T] 



point of the increase in brightness associated with a SGS 
wall in the 2D image. Uncertainties associated with this 
step may be roughly quantified by a simple error estimate 
on the she ll wi dths. This is implemented in the analy- 
sis in |4.1| and |4.2[ by varying the defined thickness of 
each shell by ± 50% of its nominal value. A related con- 
cern is shell 'leakage'. In a highly structured ISM, wind 
and supernova feedback may leak beyond the edges of a 



Chen et al. 


2000 


Walch et al. 


2012) 



However, this is 

not likely to be a major issue m the present work. We 
are concerned with identifying regions where the ISM has 
been swept-up into dense structures by large-scale stellar 
feedback. Our method of defining SGSs directly recov- 
ers this behavior, by returning the edges of dense walls, 
making no ad-hoc assumptions about the morphology of 
a shell. If there is material flowing through breaks in shell 
walls, it will either be sufficiently energetic (and/or suf- 
ficiently confined) to excavate a new dense-walled cavity, 
or it will dissipate without strongly affecting the configu- 
ration of the neutral ISM. In the former case our method 
will accurately account for it, provided the new cavity 



is well joined to the original shell. Indeed, this kind 
of phenomenon can be seen in the case of Complex 8, 
which has several lobes extending from the main body 
of the SGS. Conversely, if the outflow is unable to in- 
fluence the neutral gas sufficiently to form a void/wall 
configuration then its influence on the molecular cloud 
formation process is almost certainly negligible, and we 
are unconcerned with it in this work. We are therefore 
confidant that our shell definition method recovers the 
relevant material for our analysis. 

The presence of smaller shells within the defined SGS 
and background zones may also potentially affect the 
analysis. If such shells significantly affect the molec- 
ular gas fraction, then a systematic overabundance of 
these objects in either SGS or background zones would 
bias our result. Figure [3] overplots the Hi supershells 
of KDSB99 on a map of the SGS and local background 
zones, with Hn regions identified in the literature as su- 
perbubbles also plotted. While it is possible that some 
of these Hi objects are not genuine feedback structures, 
they nevertheless provide a rough estimate of the distri- 
bution and abundance of smaller shells in the LMC disk. 
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Figure 4. Placement of smaller LMC supershells relative to SGSs 
and backgr o und z ones. Dashed circles are the Hi 'giant shells' of 
|Kim et al.] l |1999| ). Solid circles are Ho re gions that have been 
ident ified as superbubbles in the literature (IChu et al.||1995| lOeyl 
|1996||Dunne et al.||2001| . Shaded dark grey areas show the SUs 
zones used m the analysis, with shell widths denned as in table 
Shaded medium grey areas mark the maximum extent of the 
loc al ba ckground zones used in §4.2| The global background used 
in §4, 1| corresponds to the shaded light grey and shaded medium 
grey zones. The thin dotted line marks the boundary of the region 
observed in CO. 



These shells cover very similar fractions of the total sur- 
face areas of the SGS and background zones. 0.19 ±0.01 
of the SGS zones, 0.18±0.01 of the remainder of the LMC 
disk and 0.17 ± 0.01 of the aggregate local background 
zones are occupied by small shells, where the uncertain- 
ties reflect the values over the full range of shell and back- 
ground zone widths used in the analysis. This suggests 
that to the first order, smaller shells may be assumed to 
contribute approximately equally to /h 2 in both SGS and 
background zones. Their presence is therefore unlikely to 
cause significant bias in the measurement of the effect of 
SGSs on the molecular fraction. The decision not to ex- 
plicitly include smaller shells in the analysis does impact 
on the wider astrophysical int erpr etation of our results, 
however, which we discuss in §5.1[ 

Finally, the simple method of defining local back- 
ground zones may lead to a possible systematic bias in 
the analysis of individual shells. Many SGS complexes 
arc located in the outer regions of the disk. However 
NANTEN CO coverage is limited at large galactocentric 
radii, affecting our ability to define background zones on 
all sides. As a result, background zones tend to preferen- 
tially sample material at smaller radii than their paired 
SGS zones. Since there is generally very little CO in 
the outer regions of gal axy disks (whereas H i remains 
strongly detectable, e.g. Wong 



Blitz 2002 

JL 



such an 

effect (if present) would result in values of /H 2 ,b g that 
are erroneously high, biassing the the analysis against a 
positive result. 



We sum the Hi and CO emission within all SGS zones 
to derive a global molecular fraction for supergiant shells, 
and compare this to the global molecular fraction in 
the rest of the LMC disk. We find that the total 
molecular fraction in the regions occupied by SGSs is 
/h 2 ,SGS — 0.054± 0.004. The equivalent quantity for the 
remainder of the observed area is /H 2 ,b g = 0.055 ± 0.002. 
Here, uncertainties are estimated by varying the shell 
boundaries by ±50% of the defined thicknesses, as given 
m table [TJ 

The majority of the molecular gas in the LMC is lo- 
cated in one of two large regions in which /h 2 is en- 
hanced with respect to the global average, and whose 
origins appear to be related to the global structure of 
the LMC disk. The closeness of /h 2 ,sgs and /H 2 ,b g is in 
fact somewhat serendipitous; arising from the fact that 
these regions of high molecular fraction are by coinci- 
dence equally represented in both SGS and background 
zones. The first of these CO-rich zones is the bright 
ridge-like feature sometimes known as the South-Eastern 
HI Ovcrdcnsity (SEHO), which extends for ~ 3 kpc along 
the South Eastern edge of the disk. This region is sub- 
ject to tidal forces (e.g. Staveley-Smith et al. 20031 and 
ram-pressure effe cts from the LMCs passage through the 
Milky Way Halo ( |de Boer et al.|1998l ), and hosts - 1/3 of 
the total molecular mass in the LMC. The SEHO is heav- 
ily sampled by SGS Complex 2. The second is the central 
portio n of the disk, whe re large-scale gravitation al insta- 
bility ( Yang et aL]|2007 ) and disk perturbations ( |Dottori| 



et al 
bar li 



4.1. 



4. RESULTS 
Global Molecular Gas Fraction 



996) arising from the influence of the old stellar 
cely combine to produce high gas densities. The 
arrangement of the supergiant shells is such that both 
shell and background zones sample approximately equal 
fractions of these high /h 2 regions, with the contribu- 
tion of Complex 2, which includes most of the SEHO, 
offset by the lack of SGSs in the central disk. The mean 
Hi column density is also enhanced in these two regions, 
in keeping with the results of W09, who find that the 
probability of a CO detection increases with increasing 
iVHi. _ 

This implies that SGSs are not the dominant driver of 
molecular gas production in the LMC; the global struc- 
ture of the disk is a better determinant of where molecu- 
lar matter is formed. Nevertheless, it may be that SGSs 
make some small contribution to molecular cloud forma- 
tion, that can be recovered by comparison with back- 
ground zones local to the shells. 

4.2. Molecular Gas Fractions for Individual Shell 
Complexes 

The impact of unrelated global variations in /h 2 can be 
minimized by defining background zones locally, so that 
each SGS is paired with a background zone that only 
contains material in its immediate vicinity. This kind of 
comparison has already been carri ed out successfully fo r 
selected objects in the Milky Way ( Dawson et al.|20 1 1b) 

Local background zones are defined as described in j j3.3| 
and shown in figure [3] The molecular fraction in each 
individual SGS is calculated and compared with that in 
its paired background zone. To quantify sensitivity to 
the choice of region, /h 2 ,sgs and /H 2 ,bg are computed 
for a grid of shell and background zone widths. As for 
the global case, shell boundaries are varied by ±50% of 
the defined thicknesses, as given in table [TJ Background 
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Figure 5. Box plot of log(/H 2 ,SGs//H 2 .bg) f° r individual SGS 
complexes and their local background zones. Each box illustrates 
the distribution of values obtaine d for a grid of background zones 
and shell widths, as described in §4.2| The central lines show the 
median values, the upper and lower bounds of the boxes show the 
upper and lower quartiles, and the whiskers extend out to the max- 
imum and minimum values of the data, or to 1.5 interquartile range 
if there is data beyond this range. Means are plotted as crosses 
and outliers as circles. Complex 9 contains extreme high outliers 
arising from combinations where background zones contain almost 
no detected CO; these lie outside the vertical range of this plot 
and are not marked. Each complex is color-coded by the total 
mass, Mjji + Mh 21 contained in the SGS zone. The dotted line at 
l°g(/H 2 ,SGs//H 2 bg) = marks the point at which the molecular 
fractions in the SGS and background zones are equal. The filled 
star marks the value obtained for Complex 2 (LMC2) when using 
the rest of the SEHO as the background zone (see text). 

zone widths are varied between the limits described in 



^3.3 (beginning at 10' from the shell edge and ending 
at 33' from the inner rim). Both shell and background 
zones are incremented in 2' intervals. 

Figure [5] shows a box plot of log(/ H2 ,sGs//H 2 ,b g ) for 
each SGS complex, showing the maximum and minimum 
values, the lower and upper quartiles, and the mean and 
median results obtained over the full grid of shell and 
background sizes. We find that of the nine complexes, 
five show strong evidence for increased molecular fraction 
in the SGS volumes (Complexes 2, 3, 4, 5 and 9), two 
show good evidence of a decrease (Complexes 1 and 6) 
and two show no significant difference to within the la 
uncertainties (Complexes 7 and 8). The median values 
of log(/ H2 ,sGs//H 2 ,bg) range from -0.29 to 0.57, and are 
listed in table [3j Those complexes showing an apparent 
enhancement contain ~ 70% of the total mass in the 
sample. In this context it is interesting to note that the 
complex with the strongest apparent decrease (Complex 
6) is dominated by one o f the least robust SGS detections 
in the sample (see ^3.1). 



These numbers may be used to obtain a crude esti- 
mate of the contribution made by SGSs to the formation 
of molecular gas in the LMC. We compute the difference 
between the measured Mh 2 in each SGS complex and an 
'expected' value derived from the median value of /H 2 ,bg 
- assumed to represent the molecular fraction that would 
have been measured if no SGS had been present in the 
volume. The results for each complex are given in table 
[3j In the case of Complex 2, an alternative approach is to 
use the remainder of the SEHO as the background zone, 
since there is evidence that the relationship between TVhi 
and /h 2 there may be different to the rest of the galaxy. 
This does not substantially change the results, but pro- 



vides a slightly higher estimate of /H 2 ,bg (listed in the 
table and also marked on figure [5j) . Table [2| shows the 
total Hi and H2 masses in the entire LMC aisk (within 
the CO observed region), within SGS complexes, and 
the total mass estimated to have been created due to the 
influence of SGSs. Summing over the entire sample of 
shells, this estimate suggests that ~ 25% of the molecu- 
lar matter in SGS complexes may have been formed as a 
result of the stellar feedback that formed the shells. This 
is equivalent to ~ 11% of the total molecular mass of the 
LMC. These figures fall to - 14% and - 8% when the 
background zone for Complex 2 is defined as the remain- 
der of the SEHO. Although it is the contribution from 
Complex 2 that dominates these numbers, we note that a 
positive result (~ 12% of the CO in SGSs, corresponding 
to ~ 4% of the total molecular mass in the disk) is still 
found when the SEHO is completely excluded from the 
analysis. 

We stress that these are crude estimates. Nevertheless, 
the results suggest that the formation of supergiant shells 
has a measurable positive effect on the molecular fraction 
of the LMC, albeit one that may be small compared to 
other drivers of molecular gas production in the galaxy. 
These numbers are also b est c onsidered as lower limits, 
as we will discuss below (j|5.1[). 



4.3. CO Detection Fraction in Supergiant Shells 

W09 find that the probability of detecting CO at any 
given position in the LMC disk increases with increasing 
iVjji (although it never reaches unity), and that where 
CO is present there is also a weak positive correlation 
between iVm and Nr 2 ■ Given this, the enhanced /h 2 in 
SGS zones might be expected to arise from a tendency for 
them to possess a higher fraction of pixels at the largest 
Hi column densities. However, this does not appear to be 
the case. The top panel of figure [6] shows Hi column den- 
sity histograms for the SGS and local background zones 
(aggregated). The SEHO has been excluded, since the 
relationship between CO and Hi there is different to the 
rest of the disk (W09). The shells show a low column den- 
sity tail corresponding to their evacuated voids, and their 
overall JVhi distribution is in fact skewed towards slightly 
lower values than the aggregate backgrounds. The figure 
illustrates the results for the largest background zones 
used in the analysis, but the key features of the distri- 
bution are insensitive to the choice of background zone 
width. The result is also not strongly driven by any par- 
ticular complex and holds true on individual scales for 
all SGS and background zone pairs, with the exception 
of Complex 5. 

Instead, the SGS volumes show an enhanced CO de- 
tection fraction at high Nsi- The lower panel of figure [6] 
shows the CO detection fraction as a function of Hi coh 
umn density for the same aggregate SGS and background 
zones. The CO detection threshold for both samples is 
similar, at ~ 5 x 10 20 cm~ 2 . However, the SGS zones 
show a significantly enhanced likelihood of CO detection 
for JVhi > 3 x 10 21 cm~ 2 . This behavior is mostly driven 
by several bright CO/Hi features in the dense walls of 
Complexes 3, 4 and 5 - three of the five objects show- 
ing individually enhanced molecular fractions. Similarly, 
the lower detection fractions observed at high Nm in 
the aggregate backgrounds are driven by a few patches 
of bright Hi with no CO. The small number of regions 
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Table 2 

Total Hi and H2 masses in various regions of the LMC. 



Mhi/M M„ 2 /M q 

Whole LMC, where CO observed 3.0 (2.1) b X 10 s 1.7 (l.l) b X 10 7 

SGS complexes, where CO observed 1.2 (0.8) b X 10 s 7.0 (3.3) b X 10 6 

Estimated mass created due to SGSs a ■ ■ ■ 1.8 [1.4] c (0.4) b x 10' 



a See 34. 2| 

b Values in round brackets obtained when 
c Value in square brackets obtained when 
for Complex 2. 



SEHO completely excluded from analysis, 
remainder of SEHO used as local background 
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Figure 6. Upper panel: TVjji histograms for the aggregate of SGS 
zones (light grey) and the aggregate of their local background zones 
(black line), excluding the SEHO. Dark shading and dotted lines 
show the distribution of pixels for which CO was detected in the 
SGS zones and background zones respectively. The results dis- 
played here use the maximum background zone width (see j ]4.2| , 
but the basic features of the distribution remain unchanged over 
the full range of widths tested. Lower panel: CO detection frac- 
tion histograms for aggregate of SGS zones (light grey) and the 
aggregate of their local background zones (black line), excluding 
the SEHO. Error bars are the scaled l/\/N counting errors. 

driving this trend, and the large scatter observed gen- 
erally in the CO/Hi relationships throughout the LMC, 
urges caution in interpretation of this behavior. Never- 
theless, it does suggest that the simple explanation of 
SGSs possessing higher Hi column densities and there- 
fore more molecular gas does not satisfactorily explain 
the enhanced /h 2 in supergiant shells. 

5. DISCUSSION 

5.1. The Role of Stellar Feedback in Molecular Cloud 
Formation 

It has been proposed that while the feedback from mas- 
sive stars plays a role in structuring the di sk ISM (e.g. 
de Avillez fc Berry|200l) |Joung et al.||2009| [Dobbs et al.| 



2011 Hill et al.||2012j ), and is likely responsible for trig- 



on local scales (e.g. |Boss||1995 Yamaguchi et al.||1999[ 



Hosokawa & Inutsuka 2006: Dale et al. 2007: Lcao et al 

u. , 11 , 1 1_ 



2009| Deharveng et al. 2010| ), the initial production of 



dense (molecular) clumps and clouds on galactic scales 
is driven mainly by a combinatio n of global gravitational 



instability and turbulence (e.g. [Wada fc Norman 2001 
Elmegreen|[2002l ITasker fe lhn||2TO[ |Bournaud et al' 



hlmcgrc 

2010| Mac Low" fc Glover| |2012| t Our results appear 
broadly consistent with this picture. While supergiant 
shells are clearly important in structuring the LMC disk, 
the results above suggest that they are not the dominant 
driver of the formation of molecular clouds. This is in- 
ferred both from the comparison of SGS volumes with the 
rest of the disk - which shows that the global structure 
of the disk is a better determinant of where molecular 
matter is found - and from the comparison of SGSs with 
their local environment - which finds an enhancement in 
molecular fraction corresponding to only ~ 4-11% of the 
total molecular mass of the LMC disk. 

An important caveat to the above interpretation is 
that the numbers derived in { 4.2 are (approximate) lower 
limits on the total contribution of stellar feedback to 
molecular cloud formation in the LMC. This is firstly 
due to the decision to ignore smaller supershells, which 
may also be playing a role in molecular cloud forma- 
tion, but whose contribution is not directly investigated 
in this work. Secondly, our analysis is limited to a cur- 
rent "snapshot" of the LMC. If molecular clouds were 
formed and dissipated within the lifetimes of the shells 
- as may be the case in objects formed through multi 
pie generations of propagating star formation (e.g. Do 
pita et al. 1985 ) - then such information is unrecover 



able within the present work. The lifetime of a typical 
GMC in th e LMC is roughly estimated by [Kawamura 



et al. (2009) to be 20-30 Myr. Dynamical age estimates 



place most LMC supergiant shells at ages of t < 10 Myr 



( Kim et aL~]|1999 ) , although such estimates are necessar 



ily crude. Estimates based on stellar population studies 



ically 10-20 Myr 




Dopita et al. 


1985 


Points et al. 


1999 


Glatt et al. 


2010 


). Thus the typical evolutionary 



germg a significant traction of observed star formation 



GMC, and it seems unlikely that multiple generations of 
molecular clouds will form and die within the lifetime 
of a supergiant shell. However, if indeed GMC lifetimes 
are ~ 1.5 times longer than the observable lifetime of an 
SGS, then clouds formed in shell walls will remain even 
after their parent shells have disappeared. Such clouds 
could fall in background zones, despite being feedback- 
formed. This would bias our measurement towards erro- 
neously low values for the fraction of molecular matter 
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formed by SGS activity. 

5.2. Final Comments: Dense Gas Tracers and Galaxy 

Types 

The quantity we have sought in this work is the fraction 
of the raw material for star formation that is produced 
in SGSs. At present, observational limitations mean that 
we are restricted to 12 CO as a tracer, which cannot dis- 
tinguish between diffuse molecular gas (that may be qui- 
escent), and the very dense material that is directly im- 
plicated in star formation. While some studies have ex- 



Chin et al. 


1997 


Wong et al. 


2006; 



| Wang et af?p 009), 



jects. However, future observations will be able to target 
higher density gas tracers more widely, to directly probe 
the material that i s most tightly correlate d with the star 

2012] and references 



formation rate (see Lada et al. 2010 



within). Seale et al? 12012) have recently reported spa- 



tially resolved observations of several LMC star forming 
GMCs in the high-density tracers HCN and HCO+. It 
will be interesting to see how the enhancement in /h 2 in 
supergiant shells is affected when only the densest molec- 
ular gas is considered. 

It will also be of interest to examine the relationship 
between supershells and molecular cloud formation in dif- 
ferent classes of galaxies. The LMC i s a dwarf irregular 
with low shear (Wei dner et al. 2010|), a larger Hi scale 
height than spirals ( Brinks et al. 2002 ) and only weak 
spiral structure. Low shear means shells can grow larger 
before they are deformed, high scale heights mean they 
can expand further before vertical blowout and depres- 
surization, and the lack of a strong spiral potential mit- 
igates against disruption by spiral arms. While Milky 
Way work suggests a significant quantit y of molecular 
gas for mation in two Galactic supershells ( Dawson et al. 
2011b), no systematic study of multiple shells has yet 
been performed to explore whether this behavior exists 
on large scales. Similarly, no quantitative comparison of 
CO and Hi in shells in other galaxies has been made to 
date. 

6. SUMMARY AND CONCLUSIONS 

The large-scale stellar feedback from OB clusters is 
one of several astrophysical processes postulated to drive 
the formation of molecular clouds. We have used Hi 
and 12 CO(J=1-0) data to examine the role of supergiant 
shells in the formation of molecular gas in the LMC, by 
comparing the molecular fraction, fn 2 , in SGS volumes 
to that in undisturbed 'background' regions. 

We have summed the Hi and CO emission over all SGS 
volumes to derive a global molecular fraction for super- 
giant shells, yielding /h 2 ,sgs = 0.054 ± 0.004. This is 
identical (within the errors) to the molecular fraction in 
the remainder of the LMC disk. This suggests that su- 
pergiant shells are not a dominant driver of molecular 
cloud formation in the LMC. Indeed, the global struc- 
ture of the LMC disk is generally a better determinant 
of where the highest molecular fractions are found. 

However when the impact of galaxy-scale variations in 
/h, is minimized, we find that supergiant shells have a 
positive effect on the molecular gas fraction in the vol- 
umes of space they occupy. This is measured by com- 
paring SGSs with their local surroundings, which are a 



better proxy for the undisturbed medium of individual 
objects. In this local analysis the majority of SGSs (5 
out of 9 complexes, equal to ~ 70% by mass) exhibit 
molecular fractions that are significantly enhanced with 
respect to their local background regions. Averaged over 
the full population, our results imply that ~ 12-25% 
of the molecular mass in supergiant shell systems was 
formed as a direct result of their action on the ISM. This 
corresponds to ~ 4-11% of the total molecular mass of 
the LMC. 

These figures are an approximate lower limit to the 
total contribution of stellar feedback to molecular cloud 
formation in the LMC, and constitute one of the first 
quantitative measurements of feedback-triggered molec- 
ular cloud formation in a galactic system. The results 
of this work are tentatively consistent with a scenario 
in which large-scale gravitational instabilities - as well 
as other global dynamical processes such as tidal forces 
and ram pressure in the South Eastern regions - are 
responsible for the majority of dense gas formation in 
the LMC, with stellar feedback making an important 
secondary contribution. However, a detailed analysis, 
including smaller feedback structures and considering 
quantitatively the relative observable lifetimes of molec- 
ular clouds and shells, should be undertaken before this 
conclusion can be argued strongly. 
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